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In this report, we apply Optimal Control Theory to design minimum energy 7r/2 and n pulses for 
Bloch equations, in the case where transverse relaxation rate is much larger than longitudinal so the 
later can be neglected. Using Pontryagin's Maximum Principle, we derive an optimal feedback law 
and subsequently use it to obtain analytical expressions for the energy and duration of the optimal 
pulses. 
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PACS numbers: 33.25.+k, 02.30.Yy 

Optimal Control Theory [l| has been extensively used 
recently for the design of pulses that optimize the per- 
formance of various Nuclear Magnetic Resonance (NMR) 
and Quantum Optics systems limited by the presence of 
relaxation In this report we use it to derive mini- 

mum energy it/2 and n pulses for Bloch equations, in the 
case where transverse relaxation dominates. 

The Bloch equations, in a resonant rotating frame and 
when longitudinal relaxation is neglected are fill ] 
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M x = —RM X - uj y M z 
M y = -RM y + uj x M z 

where M = (M x , M y , M z ) is the magnetization vector, 
uj x ,Lu y are the transverse components of the magnetic 
field and R > is the transverse relaxation rate. If we 
make the following change of variables, see Fig. [T] 

a = ln[M(f)/M(0)] 



tan = JM% + M$/M Z 
tan0 = M y /M x 



where M = J Ml 
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Mg, we obtain 



a = -Rsm 2 9 

9 = lu± — R sin 9 cos ( 

<f> — u)u cot 9 



(1) 
(2) 
(3) 



where lu±_ = uj x sin <fi — to y cos cj>, W|| = uj x cos cj> + uj y sin <f> 
are the components of transverse magnetic field perpen- 
dicular and parallel to = (M x ,M y ), respectively. 
Note that uiu does not affect the angle 9 of the pulse. It 
just rotates M around z-axis, resulting in a waste of en- 
ergy. Thus, optimality requires wii = =>■ 4> — constant. 
Equations ^ and ([2]) are sufficient to describe the rota- 
tion and from now on we use to to denote uij_, see Fig. [T] 
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FIG. 1: The optimal transverse magnetic field ui is perpen- 
dicular to M± and its phase is constant. Without loss of 
generality, the experimental setup can be arranged such that 
to || x-axis. In this case, <f> = 7r/2 and M rotates in t/z-plane. 



Observe that when u) is unbounded, we can rotate 9 in- 
stantaneously to the desired final value 9(T) = 7r/2or7T 
without losses in a, i.e. with a(T) = a(0) = (equiv- 
alently M(T) — M(0)). This transfer requires an infi- 
nite amount of energy so it is unrealistic. A meaning- 
ful optimization problem is the following: For a spec- 
ified final value a(T) < a(0) = (equivalently speci- 
fied MIT) < M(0)), what is the optimal control u(t), 
< t < T, that accomplishes the transfer (a(0) = 
0,9(0) = 0) -> (a(T),9(T) = vr/2or7r), while minimiz- 
ing energy J T uj 2 (t)/2dtl The control Hamiltonian Q 
for this problem is 

H = -uj 1 /2 + \ g {uj - i? sin 6> cos (9) - A a i?sin 2 9 (4) 

where Ae,A a are the Lagrange multipliers. According to 
Pontryagin's maximum principle [l[ , necessary conditions 
for optimality of (u>(t), a(t), 6(t), \$(t), A (i)) are 

&H/<&u = 0=>UJ = \g (5) 
A e = = X e Rcos29 + \ a Rsin29 (6) 

A a = -ftH/tia = => A a = constant (7) 

Additionally, the optimal (ui, a, 6, \g, \ a ) satisfies [l| 

H(u,a,9,\g,\ a ) = Q, 0<t<T. (8) 
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Using (O,© the above condition becomes 

\l - 2A e i?sin0cos0 - 2A a i?sin 2 = 



(9) 



Note that when 0(t^^ 2 ) = 7r /2 (there is such a time t_/2 
for both final values 9{T) that we consider here), then (j9]) 
gives A a = Ag(i rr / 2 )/2i? > 0. Solving the above quadratic 
equation for Xg and using ([5]), we find the optimal control 
lj. Note that only the positive solution of the quadratic 
equation has physical meaning (corresponds to increasing 
0) for the cases that we study here. The optimal oj is 
given by the following feedback law 



w(0) = i?sin0(cos0+ \/cos 2 0-r-K 2 ) 



(10) 



where k 2 = 2X a /R = \j(t 7r/2 )/R 2 - Using dTDJ), the valid- 
ity of ([6|) can be easily verified. Inserting (fT0|) in ([2]) we 
obtain the differential equation for the optimal trajectory 



6 = R sin 9 y/cos 2 6 + n 2 (11) 
Integrating (p} from t = to the final time i = T we get 



sin 2 0di = 



q(T) _ 1 M(0) 
i? i? M(T) 



(12) 



If we use (|lip to change the integration from time to 
angle, we obtain 



6(T) 



sm( 



9(0) V COS 2 + K 2 



: d0 = -lnr 



(13) 



where r = M(T)/M(0) < 1. This condition determines 
k and the results for the tt/2 and tt pulses are 
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(14) 



The duration of the optimal pulses is determined from 
the relation 



T : 



,t /.e(T) 
/ / 
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(o) 0(0) 
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(15) 



using again (|lip. Note that since = 0, n are equilibrium 
points for Ijlip. we actually start from a small positive 
initial value 0(0) = e for both cases, and additionally 
for the 7r pulse we end to the value 9(T) = it — e. The 
duration of the optimal pulses as e — > is 
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(16) 
(17) 



Finally, the energy of the optimal pulses is calculated 
from 



dt = 



8(T) 



uA9)_ 
8(0) 20(0) 



rf0 



(18) 
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(a) Optimal tt/2 puis 
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(b) Optimal n pulse 




(c) Optimal trajectory 
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FIG. 2: Minimum energy tt/2 (panel a) and tt (panel b) pulses 



for M(0) = M , M(T) = 0.6M , e = 10~ J 
trajectories are also shown (panels c,d). 



The corresponding 
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FIG. 3: The energy of the optimal pulses as a function of the 
ratio r = M(T)/M(0). 



using ([10} and ([IT]). The result is 



l-r 2 



1-r 



(19) 
and 



Observe that for r — > 1 we have T w /2,T n 
oo, as mentioned above. 
In Fig. [2] we plot the optimal 7r/2 and tt pulses for 
M(0) = M ,M(T) = 0.6M ,e = 10" 3 , as well as the 
corresponding trajectories of normalized magnetization 
vector. In Fig. [3] we plot the energy of the pulses as a 
function of the ratio r = M(T)/M(0). 

To conclude, in this report we calculated minimum en- 
ergy tt/2 and tt pulses for Bloch equations in the case 
where transverse relaxation dominates, using Optimal 
Control Theory. We expect this analytical work to serve 
as a reference for numerical studies of more complicated 
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and realistic situations that incorporate for example Ion- Acknowledgements: The author would like to thank 
gitudinal relaxation and magnetic field inhomogeneity. Jr-Shin Li for his constant support. 
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